Decoherence and disorder in quantum walks: Prom ballistic spread to localization 
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We investigate the impact of decoherence and static disorder on the dynamics of quantum particles 
moving in a periodic lattice. Our experiment relies on the photonic implementation of a one- 
dimensional quantum walk. The pure quantum evolution is characterized by a ballistic spread of a 
photon's wave packet along 28 steps. By applying controlled time-dependent operations we simulate 
three different environmental influences on the system, resulting in a fast ballistic spread, a diffusive 
classical walk and the first Anderson localization in a discrete quantum walk architecture. 

PACS numbers: 03.65.Yz, 05.40.Fb, 71.23.-k. 71.55Jv, 03.67.Ac 



Random walks describe the probabilistic evolution of a 
classical particle in a structured space resulting in a dif- 
fusive transport. In contrast, endowing the walker with 
quantum mechanical properties typically leads to a bal- 
listic spread of the particle's wave function [l|. The co- 
herent nature of quantum walks has been theoretically 
explored, providing interesting results for a wide range 
of applications. They state not only a universal platform 
for quantum computing [2| but also constitute a powerful 
tool for modelling biological systems (H-IU, thus hinting 
towards the mechanism of energy transfer in photosyn- 
thesis. Quantum walks of single particles on a line have 
been experimentally realized in several systems, e.g. with 
trapped atoms [f| and ions 0, Q ; energy levels in NMR 
schemes [3, [lfl] ; photons in waveguide structures , a 
beam splitter array 12), and in a fiber loop configura- 
tion [l3| . Although these experiments opened up a new 
route to higher dimensional quantum systems, more so- 
phisticated quantum walks need to be implemented to 
pursue the realm of real applications. A first step in this 
direction has been recently reported in which two 
particles execute a simultaneous walk and display intrin- 
sic quantum correlations. 

One of the most important requirements for realizing 
quantum walk-based protocols is the ability to control the 
dynamics of the walk, that is to access and manipulate 
the walker's state in a position dependent way [ljl [l(|. In 
this paper we present the first experimental realization of 
quantum walks with tunable dynamics. We investigate 
the evolution of quantum particles moving in a discrete 
environment presenting static and dynamic disorders. 

As predicted by Anderson in 1958 [l7|, static disorder 
leads to an absence of diffusion and the wave function 
of the particle becomes localized, which, e.g. would ren- 
der a conductor to behave as an insulator. Anderson 
localization has been experimentally investigated in dif- 
ferent physical scenarios, e.g. employing photons mov- 
ing in semiconductor powders [HI and photonic lattices 
19, 20], or even via Bosc-Einstein condensates 21, 221. 



However, although theoretically predicted in the context 
of quantum walks (23l425| . the effect has never been ob- 
served in a discrete quantum walk scenario. 

Furthermore, it is interesting to note that the energy 
transport in photosynthetic light-harvesting systems is 
influenced by both, static and dynamic disorders, and 
it is precisely the interplay between the two effects that 
lead to the highly efficient transfer in those molecular 
complexes 0, Sj- Thus, in order to simulate a realistic 
influence of the environment, we go further in our stud- 
ies by investigating the effect of dynamical noise, which 
typically induces decoherence [2(3, l27[ . Utilizing the abil- 
ity to easily tune the conditions for the quantum walk, 
we demonstrate here the diverse dynamics of quantum 
particles propagating in these different systems. 

In our experiment we realize the quantum walk of pho- 
tons by employing a linear optical network. The evolu- 
tion of the particle's wave function \tp(x)) is given by 

\i;(x)}^ 7x \i>(x))+J2P*,Mk)), (1) 

with the position dependent amplitudes 7^ and /3 X (C de- 
termining the probability of the particle to stay at the 
discrete position x or evolve to the adjacent sites k, re- 
spectively. 

We study the expansion of the particle's wave packet in 
four different scenarios, (i) First of all we implement the 
quantum walk in a homogeneous lattice, showing that it 
presents an evolution that is free from decoherence. (ii) 
Next, we introduce static disorder by manipulating the 
lattice parameters 7 X and f3 Xik , thus observing Anderson 
localization. We then examine two scenarios leading to 
decoherence, which essentially differ in the time scales of 
the occurring dynamic perturbations, (iii) In this case a 
dynamic randomization of the lattice parameters 7 X and 
/3 x>h simulates the evolution of a particle interacting with 
a fast fluctuating environment. The resulting dephasing 
suppresses the underlying interference effects and hence 
causes the particle to evolve just like in a classical ran- 
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dom walk [6j, |8j. (iv) In the last scenario we simulate 
a slowly changing homogeneous environment. While j m 
and fi Xtk are stable during a single realization, a slow drift 
leads to different conditions for subsequent particles, thus 
affecting results obtained in an ensemble measurement. 

In a discrete quantum walk the position of a parti- 
cle evolves according to its internal coin state c). For 
our photonic implementation we use the linear horizon- 
tal | if) = (1,0) T and vertical \V) = (0, 1) T polarization 
of light. The state of the photon after N steps of the 
walk is found by applying the unitary transformation 
U = nr=i SC* to the initial state \tp(x) ) — \x ) ® \c ). 
The coin operation C n (x) manipulates the polarization 
of the photon in dependence on the position x and the 
step number n. In the basis {\H), \ V)} the coin operator 
is given in matrix form by 

CM - ( el * H(X> ^ f cos ( 2 ^) sin ( 20 ) ^ 

e**vW J ^sin(20) - cos(20) )' [ > 

with the diagonal matrix representing a phase shift (j) H (x) 
for horizontal and 4>v(x) for vertical polarizations, while 
the second matrix corresponds to a polarization rotation 
of 26. The step operation S shifts the position x of the 
photon by +1 if the polarization is horizontal and by — 1 
if it is vertical. 

Following Eq.([T]), the evolution of the wave function 
with the step number n is given by 

|V(x)» +a >=7-l^(a5)»)+j8-.-± a (|V(K + 2)n) + |^(«-2) B »- 

(3) 

Note that the transition coefficients 7 X and /3 X , X±2 are 
fully set by the coin operations C n+1 (x) and C n+2 (x). By 
changing the parameters <fi(x) H/v and 6 in a controlled 
way we can alter the coefficients and hence create diverse 
types of physical conditions for a quantum walk scenario. 

A simple measure to quantify the spread of the wave 
function in the different systems is provided by the vari- 
ance a 2 of the final spatial distribution. While the de- 
coherence free quantum walk presents a ballistic spread, 
with a 2 oc n 2 , the classical random walk is diffusive, char- 
acterized by a 2 — n. In contrast to both, in a one dimen- 
sional system with static disorder the wave packet shows 
exponential localization after a short initial expansion. 
The stagnation of the wave packet spread is thus evi- 
denced by a constant variance. 

The functional principle of our experimental setup is 
sketched in Fig. QTa) and is discussed in detail in jliij ]. We 
generate the input photons with a pulsed diode laser with 
a central wavelength of 805 nm, a pulse width of 88 ps and 
a repetition rate of 110 kHz. The initial polarization state 
of the photons is prepared with retardation plates. Each 
coin operation consists of a polarization rotation, which is 
realized with a half- wave plate (HWP), and a subsequent 
phase shift implemented by a fast switching electro-optic 
modulator (EOM), as described in Eq.©. The proper- 
ties of the EOM impose that (j) v (x)/(f> H (x) ~ 3.5. The 
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FIG. 1: (a) Schematic setup. The coin state is manipu- 
lated via half-wave plate (HWP) and electro-optic modula- 
tor (EOM). The two polarizing beam splitters (PBS) allow 
to implement the step operation. APDs: Avalanche photodi- 
odes; (b) Probability distribution after 28 steps of a symmet- 
ric Hadamard walk with initial circular polarization. Stacked 
bars: Adapted theory splitted into the two coin states \V) 
(blue, bottom) and \H) (red, top). Gray dots show experi- 
mental data for vertical polarization, black dots the sum of 
both polarizations. Error bars correspond to statistical errors. 



step operation is realized in the time domain via two po- 
larizing beam splitters (PBS) and a fiber delay line, in 
which horizontally polarized light follows a longer path 
(Fig. QJa)). The resulting temporal difference of 5.9 ns 
between both polarization components corresponds to a 
step in the spatial domain of x ± 1. After a full evolu- 
tion the photon wave packet is distributed over several 
discrete spatial positions or, equivalently, over respective 
time windows. For detection the photon gets coupled out 
of the loop by a beam splitter with a probability of 12% 
per step. We employ two avalanche photodiodes (APD) 
to measure the photon's time and polarization properties, 
which gives information about the number of steps, the 
specific position of the photon, as well as its coin state. 
The probability that a photon undergoes a full roundtrip 
without getting lost or detected is given by rj BStup = 0.55 
(0.22) without (with) the EOM and the detection effi- 
ciency is 77 dct = 0.06 per step. To determine the statistical 
distribution of one specific step we detected more than 
10 4 events in an overall measurement time of maximally 1 
h, limited by the setup stability. This guaranteed an ab- 
solute statistical error of the assessed probability at each 
position of less than 0.01. An average photon number 
per pulse at the detected step of less than (n) < 0.003 
ensured a negligible probability of multi-photon events 
P{n > 1)/P(n = 1) < 0.02. 

(i) Homogeneous lattice. — In the first of our four quan- 
tum walk scenarios we investigate a homogeneous envi- 
ronment, thus testing the intrinsic coherence properties 
of the setup. The spatial distribution after 28 steps can 
be seen in Fig. [IJb). We used the initial state \ip ) = 
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|0> ® 3s(\ H ) + i \V)) and the Hadamard coin (9 = tt/8) 
at each position. The final state clearly shows the charac- 
teristic shape of a fully coherent quantum walk: the two 
pronounced side peaks and the low probability around 
the initial position. Moreover, the polarization analysis 
confirms the expected dependence of the particle's final 
position on its coin state. An adapted theory includ- 
ing only small imperfections of the coin parameter 9, the 
initial coin state and differential losses between the two 
polarizations fully explains the final spatial and polariza- 
tion distribution. The quality of the result can be quan- 
tified by the distance d(P m , P th ) = \ J2 X \P m (x) ~ P th {x)\ 
between the measured P m and the theoretical P th prob- 
ability distributions. It ranges between for identical 
distributions and 1 for a complete mismatch. The dis- 
tance of the measured walk to the adapted quantum the- 
ory is d(P m ,P qw ) = 0.052 ± 0.015. For comparison we 
calculated the distance to the fully decoherent (classical) 
scenario, obtaining d(P m ,P cl ) = 0.661 ± 0.015. Hence, 
our result confirms an almost decoherence free evolution 
after 28 steps. 

(ii) Static disorder. — We implemented the evolution 
of a particle in an environment with static disorder us- 
ing a quantum walk with variable coin operation. To 
create a static disorder a coin operation is required, 
which is position and not step dependent. In our sys- 
tem this is realized by a controlled phase shift 4> H/V (x), 
such that the photon acquires the same phase any in- 
stance it appears at position x. To generate a random 
static phase pattern we applied a periodic noise signal 
to the EOM. The periodicity of the signal was carefully 
adjusted to ensure that the applied phase shift opera- 
tion is strictly position dependent. Using different phase 
patterns at subsequent runs allows to average over var- 
ious disorders, as considered in the model of Anderson. 
The strength of disorder is determined by the maximal 
applied phase shift 3> m<Me , which defines the uniform in- 
terval (f> v (x) G [-$„,„,$„,„], from which the phases 
are chosen. The probability distribution after eleven 
steps is shown in Fig. [2ja). We used the initial state 
V-'o) = |0) ® \H), 9 = tt/8 and a high disorder strength 
(f mlI = (1.14 ± 0.05)7r). In contrast to the decoherence 
free quantum walk ($ mal = 0, inset of Fig. GJc)), in the 
disordered scenario the expansion of the wave packet is 
highly suppressed. We observe a strictly enhanced ar- 
rival probability around the initial position, which also 
displays the predicted exponential decay. This striking 
signature of Anderson localization is emphasized by lin- 
ear fits in the semilog scaled plot (inset of Fig. [H[a)). 
Our results are in agreement with a theoretical model 
determined by a Monte Carlo simulation of 10 4 different 
phase patterns compatible with our experiment. Com- 
pared to (i), the number of steps is reduced due to the 
additional losses introduced by the EOM. 

(iii) Fast fluctuations. — To generate a system with dy- 
namic disorder we detuned the temporal length of the 
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FIG. 2: Measured probability distribution (front) and respec- 
tive theory (back, gray bars) of 11 steps of a quantum walk 
(9 = 7r/8) with static disorder (a), dynamic disorder (b) and 
in a decoherence free environment (Inset (c)). The insets in 
(a) and (b) show the measured distribution in semilog scale 
with linear (a) and parabolic fit (b). (c) Transition of the 
variance from ballistic quantum walk to diffusive/ localized 
evolution due to dynamic (red squares) and static (green dots) 
disorder with increasing disorder strength $ raal ; dashed lines: 
theory without adaption for experimental imperfections. The 
red solid line marks the variance of a classical random walk. 
(Vertical error is smaller than the dotsize). (d) Relative fre- 
quency /( 4>v |) of the applied phases 4>v for the signal with 
interval Q m ax = (1.02±0.05)7r. The dashed line indicates the 
uniform distribution. 



noise signal, thus eliminating position dependent phase 
correlations. Decoherence appears as a consequence of 
the dynamically varying phase suffered by the quantum 
particle during the evolution. As a result, the photon 
undergoes a classical random walk, revealing a binomial 
probability distribution (Fig. [UJb)). In contrast to the 
previous case, the spatial profile of the wave packet shows 
a parabolic shape in the semilog scale (inset, Fig. [2][b)). 

A stepwise increase of the disorder strength $ mo[c nicely 
demonstrates the controlled transition of the system from 
the ballistic evolution (decoherence free quantum walk) 
towards the diffusive evolution/localization in a scenario 
with dynamic/static disorder (Fig. Htc)). For this pur- 
pose we characterize the resulting expansion profile by 
its variance a 2 . Without decoherence (& ma:c = 0) the 
ballistically spreading wave packet shows a large expan- 
sion induced by quantum interference after eleven steps. 
The slightly higher variance in comparison to the theory 
can be explained by small polarization rotations in the 
EOM. In a system with dynamic disorder, decoherence 
reduces the expansion of the wave packet to the level of 
a diffusive classical particle. In contrast, static disorder 
leads to a stagnation of the spread and hence an even 
smaller variance. Our results clearly demonstrate how 
the amount and kind of disorder influence the expansion 
of the particle's wave packet. 
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FIG. 3: (a) Averaged probability distribution in a slow deco- 
herence scenario with different coin angles 9 € [0, 7r/4]: Mea- 
surement (orange, front) and theory (gray, back); (b)+(c) 
Trend of the variance with the number of steps: Measurement 
up to 12 steps (b) and simulation up to 50 steps (c). While 
photons undergoing the decoherence free quantum walk (blue 
triangles) and the evolution with slow dynamic disorder (or- 
ange dots) show a ballistic behavior, in a classical random 
walk (red squares) they move diffusively, and, finally, in the 
case of static disorder (green diamonds) they stagnate. 



The agreement between theory and measurement in 
the completely dephased scenario (Fig. HJb)) confirms a 
sufficient randomness of the applied noise signal. Fur- 
thermore, an independent interferometric measurement 
revealed the relative frequency of the used phases f(\(j> v \), 
as can be seen in Fig. [Sfd) with $ max = (1.02 ± 0.05)tt. 
However, small imperfections of the EOM lead to a dis- 
tribution slightly off uniformity, which also influences the 
measured variances shown in Fig. [U[c). 

(iv) Slow fluctuations. — As the fourth scenario we sim- 
ulated fluctuations in a homogeneous system, but with 
parameters that change in a time scale much larger than 
the full duration of a single quantum walk. Although the 
individual evolution is not affected under these circum- 
stances, an ensemble measurement of subsequent walks 
results in an average over coherent evolutions in different 
types of lattices. For this purpose we changed the pa- 
rameter 8 £ [0, tt/4] in steps of 7r/18 for a quantum walk 
with initial state |-0o) = |0) ® -+=(\H) + i\V)). An average 
over the full range 9 G [0,7r/4] exhibits a nearly uniform 
spatial distribution of the wave packet with an enhanced 
probability to arrive at its initial position x = after 10 
steps (Fig. EJa)). Especially the high chance to reach the 
outermost positions x — ±10 differs significantly from all 
previous scenarios. This increases the variance of the 
distribution (a 2 w) = 40.00 ± 0.42) to a level, which is 
even higher than in the decoherence free quantum walk 
with the Hadamard coin (cr 2 {i) = 31.27±0.19). The result 
demonstrates that special kinds of decoherences can even 
speed up the expansion of wave packets in homogeneous 
lattices. 



Finally, the geometry of the setup allows to observe 
easily the wave packet's evolution step by step in all four 
scenarios (Fig. EJb)). For cases (i) and (iv) we observe a 
ballistic spread, with an even faster expansion in a system 
with slow fluctuations. The evolution with fast dynamic 
disorder (iii) is clearly diffusive. Lastly, under the con- 
dition of static disorder (ii) the variance saturates after 
few steps and the dynamics is dominated by the effect of 
Anderson localization. For comparison, we show in Fig. 
[3jc) a theoretical plot for the evolution of the variance 
over fifty steps. The parameters used in simulation and 
experiment are equivalent to the experimental settings 
used for Figs. QJb), [l^a-b) andELa). 

In conclusion, we presented how disorder and fluctua- 
tions in a periodic lattice can influence the evolution of a 
traversing particle. We observed a fast ballistic spread for 
slowly changing lattice parameters, a diffusive spread in 
the case of dynamical disorder and Anderson localization 
for lattices with static disorder. Furthermore, we showed 
the controlled transition between the different regimes. 
The high flexibility and control allows not only the study 
of further decoherence phenomena in quantum walks but 
also to simulate specific physical scenarios of interest for 
the solid state and biophysics community. Moreover, the 
possibility to manipulate quantum walks with time de- 
pendent coin operations is a fundamental step towards 
the realization of quantum walk-based protocols. 
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